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Abstract 

Determining the physical Hilbert space is often considered the most difficuh but 
crucial part of completing the quantization of a constrained system. In such a situ- 
ation it can be more economical to use effective constraint methods, which are ex- 
tended here to relativistic systems as they arise for instance in quantum cosmology. 
By side-stepping explicit constructions of states, such tools allow one to arrive much 
more feasibly at results for physical observables at least in semiclassical regimes. Sev- 
eral questions discussed recently regarding effective equations and state properties in 
quantum cosmology, including the spreading of states and quantum back-reaction, 
are addressed by the examples studied here. 
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^ : 1 Introduction 

^ ■ One of the key issues in quantizations of fundamental theories, which due to their covariance 
0\ '. properties are systems with gauge freedom generated by constraints, is the determination 
^ ' of physical observables. They must satisfy the constraint equations and be invariant under 
gauge transformations. For canonical quantum theories, solving constraints is traditionally 
^ ' done at the state level: one constructs a physical Hilbert space of states annihilated by the 
■ constraint operator(s) and equipped with an invariant inner product. Explicit constructions 
can be done in some special cases by different methods. 

Since explicit derivations are possible only in specific cases, it is not always clear whether 
the results are generic or mere artefacts of the simple models used. It is therefore impor- 
tant to have approximate methods for a wider range of cases, or at least to be able to 
perturb around known solvable ones while still ensuring that the constraints are solved 
and the observables are gauge invariant. It turns out that such perturbation schemes are 
most feasible if one deals with the observables directly, such as expectation values, side- 
stepping the computation and physical normalization of states. This procedure gives rise 
to canonical effective equations and constraints [HE]- 

A procedure for effective constraints has been formulated in [2] and applied to param- 
eterized non-relativistic systems with a constraint pt + H = where pt is the momentum 
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of time and H the Hamiltonian. (The concepts and results are reviewed briefly below.) It 
was shown that the physical observables in suitable regimes, including semiclassical ones, 
can be determined without making use of a physical inner product but instead through 
implementing reality conditions for quantum variables such as fluctuations, correlations 
and higher moments. For applications of these methods to quantum gravity and cos- 
mology one has to extend them to relativistic systems, offering one additional subtlety: 
constraints would now be of the form — = 0, requiring one to take a square root 
and to make sign choices. Mathematically, for instance, the question is how to precisely 
define y/H'^ = \H\ at the operator level. This may not be obvious if the operator H has a 
complicated spectrum or is not positive definite. Physically, one must decide how to treat 
and separate positive and negative frequency solutions corresponding to the two solutions 
Pt = ±\H\. (See e.g. [21 H] for discussions of relativistic systems.) 

For a time-independent Hamiltonian, it turns out that one can, at least for semiclassical 
purposes, simply use pt = ±(H) as the effective Hamiltonian without absolute values, 
even if H is not positive definite. One only has to ensure that the initial values used in 
the effective equations of motion correspond to an initial state supported on a part of the 
spectrum of H with a definite sign. On such a state and with a self-adjoint H, (|H|) = (H) 
if the state is supported only on the positive part of the spectrum of H, and (|H|) = —(H) 
if the state is supported only on the negative part. Since the Hamiltonian is preserved, 
these statements will hold true throughout the whole evolution and there is no need for 
an absolute value in the effective Hamiltonian. This fact has been made use of in several 
recent derivations of effective equations in quantum cosmology, where the relevant versions 
of H are not positive definite [TJ [8] . 

In those models, deparameterization was performed using t = as an internal time 
from a free, massless scalar 0. The same types of models also allow the construction of a 
physical Hilbert space at the state level P [TU] , with the results in agreement with those 
obtained from effective equations. Most interesting from the cosmological perspective is, 
however, a system where the scalar (p has a non-trivial potential or at least a mass term. 
This has two immediate implications: in general, one can no longer deparameterize globally 
since the solutions for the scalar would not be monotonic in the time coordinates, and the 
Hamiltonian would no longer be (internal) time independent. Positivity can no longer be 
ensured just by an initial condition, and using p^ = ±(H(0)) as an effective Hamiltonian 
without an absolute value may then seem questionable. Explicitly using an absolute value, 
on the other hand, would make a derivation of the effective equations more complicated. 
At this stage, a direct treatment of effective constraints for relativistic systems without 
deparameterization becomes relevant. This is what we present in the current paper. 

We will consider in detail models of relativistic particles and properties of observables 
in their quantum theories. In the massive case, for instance, we are dealing with the 
quantization and implementation of a constraint C = p'f — p^ — rn?- For physical states, 
the expectation value 

{C)=pl + {Aptf-p'-{Apf-m'' 
must vanish and thus represents a quantum constraint. (We will notationally identify 
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classical degrees of freedom with the expectation values to simplify the notation and to 
show the relation between classical terms and quantum corrections. Thus, E = (E) and 
(A_E)^ = ((E— (E))^) = (E^)— E^.) As we will see in more detail below, there are additional 
independent constraints since expressions such as (qC) must also vanish for physical states. 
Allowing all possible factors to the left of C, this presents a constrained system of infinitely 
many constraints for infinitely many quantum variables given by the moments of a state. 
The combined system of all constraints must be solved to find observable results, which is 
feasible in semiclassical regimes where only a finite set of moments suffices to characterize a 
state approximately. The same kind of approximation also allows one to include potentials 
within the constraint, which may be explicitly time-dependent. We will exploit this to 
justify the procedures used in quantum cosmology for deparameterized systems with time- 
dependent potentials, as developed in pTl [T2| [13]. 

Another question of interest is that of the spreading of states and quantum back- 
reaction of fluctuations and higher moments on the expectation values. If we compare the 
effective constraint (C) written above, which contains only the second order moments in 
addition to the expectation values, with the effective Hamiltonian of the corresponding 
deparameterized system. 



(H) = (v^^^+^) = (V(p+(p-j9))2 + m2= V^^T7;^+^- ((p-pD 



n=2 

,2/'A^^2 ^2 



- VP +^ + ^^2 + ^2)3/2 (p2 + ^2)7/2 UP ) + 

with a whole formal Taylor series that includes higher moments, different coupling terms 
between the expectation values and the moments seem to be implied. Then, back-reaction 
might seem different in these two treatments, apparently making them incompatible. By 
our specific constructions in this paper we will reconcile these apparent disagreements, 
and provide an illustration by numerical solutions in a specific example (App. [B]). This 
is also important for quantum cosmology, where quantum back-reaction is crucial to the 
understanding of how a quantum state evolves toward and possibly through the big bang 
and how much of the pre-big bang state can be reconstructed [15] . 



2 Effective constraints 

In a canonical effective description, the dynamics of a quantum system with n degrees of 
freedom is formulated in terms of the expectation values, i.e. the evaluations of a state 
functional in the elements of an algebra generated by In basic operators and pi, 
2 = 1, . . . , n. This whole set of infinitely many variables can be conveniently split into the 
2n expectation values of the basic operators, such as gj = (qj) and = (pj), i = 1, . . . , n, 
together with the infinitely many moments of the form 

= /j](q, - (q,))»»(p, - (p,))M 

\«=1 / Weyl 
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where the subscript "Weyl" denotes the totally symmetric ordering of all factors involved. 
Between all these variables a Poisson bracket is defined following from the algebra of 
commutation relations: 

{(A),(B)} = ii^ (2) 

extended using linearity and the Leibniz rule. 

At second order, J2i=ii^i + ~ 2, this set of moments includes all fluctuations and 
covariances. In semiclassical regimes, moments fall off at least as /i2 2]i=i(ai+''«) such that 
only low orders need to be considered for the first approximation to quantum effects. 
Below we will employ the notation (Aa)^ = ((a — a)^) for fluctuations and A{ab) = 
((a — a)(b — 6))wcyi = |((a — a)(b — 6) + (b — 6)(a — a)) for covariances, which may make 
second order equations easier to interpret. 

Any operator C = C(qj,pj) gives rise to a function on the space of states, which can 
be expressed in terms of the moments by Taylor expanding 

(C(qi, pj)) = {C{q, + (qi - qi),Pj + {pj - pj))) (3) 

in qj — qi and pj — pj as in ([T]). If C(qi, pj) is a constraint operator, ([3]) must vanish on 
physical states and thus is a constraint on the quantum phase space. 

A single constraint on the phase space removes one pair of variables, but not the whole 
tower of moments associated with it in the quantum phase space. For a complete reduction 
one has to make use of additional constraints, provided by the set of phase space functions 
(/(qi, Pi)C(qj-, pa;)) which must also vanish on physical states. These functions are in 
general independent from the quantum constraint (C(qj,Pj)) as functions of expectation 
values and moments. As shown in [2j, this set of constraints remains first class. (Note that 
we do not order the products of the operators in the constraints symmetrically to ensure 
that the constraint operator acts directly on the state. As a result, in some cases one has to 
deal with complex-valued constraints requiring reality conditions for physical observables. 
This has been discussed in [2] and will also be seen in more detail in the examples below.) 

Given such a system of constraints on the quantum phase space, one can proceed in 
the classical way and find the reduced quantum phase space of observables or solve the 
constraints and fix the gauge. At this stage, it is convenient (but not required) to decide 
which internal time variables {t,pt) among the {qi,Pj) should be used. Since a quantization 
of the corresponding deparameterized system, if it exists, would not give rise to any moment 
involving an operator of t or pt, solving the constraints must remove all moments including 
at least one factor of t or p^ from the original quantum phase space. That this indeed 
happens was verified to second order of the parameterized non-relativistic particle in |2]. 
In that case, there was a single linear term pt in the constraint, such that all pf-moments 
{AptY, A{ptt), A{ptq) and A{ptp) can be removed by solving the system of constraints, to 
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second order: 

((q-g)C) = A(gp,) + |:| + |:AM = 

((t-t)C) = |:A(pt) + A(tp,) + | = 
{{Pt-Pt)C) = {Ap,Y + ^A{pp,) = 
((p-j9)C) = A{pp,) + ^{Apy = 0. 

This leaves the moments (At)^, A{tq) and (Atp), which are removed by factoring out 
the gauge flow, or simply by setting them to zero as a well-defined set of gauge-fixing 
conditions. (Note that a smaller number of gauge fixing conditions than constraints is 
required because the second order moments satisfy a Poisson algebra which is degenerate 
from the symplectic point of view; see [16] for some notions of constrained systems in the 
non-symplectic case. Additionally, setting the fluctuation (At)^ to zero is consistent with 
the generalized uncertainty relation 

(AtnAp^y - A{tp,y > ^ 

since t — prcorrelations A(tpt) = —ih/2 are required by the constraints, especially ((t — 
t)C) = with the gauge fixing condition A{pt) = 0, to be imaginary and of just the 
right size to saturate the uncertainty relation.) After solving the constraints and fixing 
the gauge, observable moments are recovered on which physical reality conditions can be 
imposed. 



3 Free relativistic particle 

Classically a free relativistic particle is described by a single constraint 

2 2 2 

C = Pf — p — m 

on the phase space coordinatized by two canonical pairs t,pt and For the quantum 

version we consider the unital associative algebra generated by four basic elements t, pt, q, p 
subject to the canonical commutator relations. That is, consists of (countable) sums of 
polynomials of the form t^p^q'^p'; terms with a different ordering may be expressed using 
the commutation relations 

[t, pt] = iht , [q, p] = ihl . 

There is no product-ordering ambiguity in the case of the above constraint and it is nat- 
urally identified with an element C = p^ — p^ — m^l of =2/. As the equivalent of Dirac's 

•^We assume the units have been chosen so that both length and momentum have the units of the square 
root of action (e.g. geometrized units). 
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condition Cijj = we demand that the constraint has a vanishing right action on the 
states, which in our case are complex hnear functions a : ^ <C, this imphes that 
a(aC) = 0, Va e henceforth we drop exphcit reference to a and write this condition 
as (aC) = for the expectation value in a physical state a. In order to impose all these 
conditions systematically we take the previously mentioned basis in ^ and impose the 
constraint via an infinite (but countable) set of conditions 

(tV^q'^p'C) = . (4) 

We reduce the above infinite system of equations using the same method that was 
previously employed for a Newtonian particle [2] — a semiclassical expansion based on the 
hierarchy 

((t - imp, - (p,)y (q - (q))^(p - (p))')^^^, (X 

of moments. 



3.1 Constraints at second order 

3 

In what follows, we assume a semiclassical state and drop the terms of order and above, 
keeping the terms of order h and below. This will suffice to demonstrate the feasibility 
of our methods for relativistic systems. To this order our system is described by fourteen 
independent functions: four expectation values of the form a = (a); four spreads of the 
form (Aa)^ = ((a — a)^) and six covariances of the form A{ab) = ((a — a)(b — 6))weyi- 
(Poisson brackets between all these variables are listed in App. \M) The infinite system of 
constraint functions is reduced to just five non-trivial conditions 

C = {C)=p^,-p^-m^ + {Aptf - {ApY = 

Ct = {{t - {t))C) = 2ptA{tpt) + ihpt - 2pA{tp) = 

Cp, = {{pt-{pt))C)=2pt{Apty-2pA{ptp)=0 

Cg = {{q-{q))C) = 2ptA{ptq)-2pA{qp)-ihp = 

Cp = {{p-{p))C) = 2ptA{ptp)-2piApf = 0. (5) 

Note that the semiclassical hierarchy of variables is critical to the above reduction in 
the number of constraint conditions. In particular, C = implies that pi — p^ — to^ = 
{ApY — (Apt)"^, which in turn implies that the combination of the expectation values 
Cciass '■= vi ~ ~ is of order h on the constraint surface. In other words, the classical 
constraint is satisfied to order h. The terms of the form (Aa)^Cciass and A(a6)Cciass are 
then of order fi^ and should be dropped in our present treatment. The complete infinite 
system of constraint functions is a closed Poisson algebra — or, in the language of classical 
constraint analysis, a first-class system [17]. In general, due to the nature of the above 
truncation one would expect the reduced system of constraints to remain closed only to 
order h. In our case the Poisson algebra of the truncated set of constraint functions — 
displayed in Table [1] — is exactly closed with respect to the bracket. 
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Table 1: Poisson algebra of constraints for a free particle. First terms in the bracket are 
labeled by rows, second terms are labeled by columns. 
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To solve the constraint system we eliminate five variables using the five conditions 
from (IS]). Specifically, we eliminate the five quantum variables associated with pj, having 
in mind that t will be chosen as time in a deparameterized treatment. We start by noting 
that ptCp^ = gives 

= Pti^ptf - pptA{ptp) . (6) 

However Cp = implies 

PtA{ptp) = p{Apf 

and substituted in ([6]) gives 

= p',iAptf-p\Ap)\ 
Finally, eliminating pf through C = we obtain a quadratic equation in {Apt 

{{Aptfy - (Api)'(/ + + {Ap)')+p\Ap)^ = 

with two solutions 

[AptY = ^(^p^ + m^ + {Apf ± ^(p2 + ^2 _^ {Apff - Ap^{Apy 

In order to see whether either solution is compatible with the hierarchy assumed by the 
semiclassical approximation, we expand the solution to order %. One finds 



2 \ p^ + m? V [p^ + m?Y 

= lip' + m^) (l ± 1 + J^^^2iP' + ± -P')) + O ((Ap)^)) . 
Looking at the solution with the "+" sign we see the following leading order behavior 

{Aptf = p^ + + 0{h) 

which is inconsistent with the assumption that {AptY is of order h. The solution with the 
"— " sign leads to 
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which is of order h and therefore consistent with the semiclassical approximation. Sub- 
stituting the latter solution back into the constraint conditions ([5]) we obtain two sets of 
solutions 



Pt = ±E 
) 

A{ptq) = ±|(^A(gp) + y 

AM = ±|(Ap)2 (7) 

where ^ 
E = ^ (^p^ + m'^ + [Apf + v/(p2 + m2 + [Ap^y - Ap^{Apy^ ' . 

By rearranging the terms in the above expression one can easily verify that E >Q when 
reality and positivity of the physical variables are imposed (see Section [32]). One recovers 
the usual relativistic dispersion relation where energy equals \/p^ "irfi if one assumes a 
"momentum eigenstate", that is if one sets the spread in momentum (Ap)^ = 0. 



3.2 Gauge freedom 

The truncated system of constraints ^ is equivalent (when consistency with the semi- 
classical approximation is evoked) to the restriction to two disjoint surfaces, each one 
corresponding to a choice of sign in ([7]). Each surface is described by an equivalent set of 
"linearized" constraints 

Ci± = pt±E 

C2± = {Aptf-p'-m^-{Apf + E^ 
C3± = A{p,p)±^{ApY 

C4± = A(p,g) ± I (A(gp) + I) 

a± = A(tp,)±|A(^P) + y- (8) 

The above constraint conditions can be expressed as sums of the conditions in ([5]) and 
therefore form a first class system. Additionally, for the calculations to follow it is useful 
to note that pf, p, {AptY, A{ptp) and (Ap)^ are first class functions with respect to either 
set of constraints. This can ultimately be traced back to the fact that [pt, C] = = [p, C]. 
It follows that E, which is a function of p and {Ap)^ only, is also first class. 

On the constraint surfaces the "linearized" constraints can be used to eliminate the 
five variables pt, {AptY, A(tpt), A{ptp), A{ptq). At this stage there remain four degrees 
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of freedom associated with the algebra elements generated by t. These will be treated 
as gauge parameters associated with the time evolution of the system. From this point 
of view, one has a four-parameter space to choose from when it comes to the evolution 
of the "physical variables" (i.e. those associated with the algebra generated by q and p 
alone). Viewing our system expanded to second order in quantum variables as a classical 
constrained system, the evolution on the "physical variables" — q, p, (Ag)^, A{qp), (Ap)^ 
— may be generated by any constraint function of the form 

Cnam = ^^/ijCj± (9) 

i 

where the multipliers fii are arbitrary functions of the physical variables. The presence of 
several constraints, all associated with the classical Hamiltonian, means that a priori there 
is no unique time parameter. Depending on the choice of gauge, any combination of t with 
moments involving t can play the role of time. 

At this stage we would like to restrict the gauge freedom down to a single parameter 
and to interpret the first class flow in the direction of t as the dynamical evolution of the 
system. This may be accomplished by introducing three gauge choices 

01 = (At)2-/i(g,p,(Ag)2,A(gp),(Ap)2)=0 

02 = A{tq)-f2{q,p,{Aqf,A{qp),{Apf) = 

03 = A{tp)~h{q,p,{Aqf,A{qp),{Apf) = 

with functions /i, /2 and /s to be determined. We define Cnam as the first class function 
([9]) that remains after the gauge conditions have been introduced. It must therefore be a 
combination of the original constraint functions as in equation (Q that in addition satisfies 

{CHam,0*}~O, 2 = 1,2,3 

where the symbol '~' denotes equality on the surface defined by imposing both the con- 
straints and the gauge conditions. A simple set of such conditions that was also used in 
Ref. [2] to recover the deparameterized dynamics of a Newtonian particle is provided by 

01 = (At)2 = 

02 = A(tg) = 

03 = A(tp) = 0. (10) 
(Again, (At)^ = is consistent with the uncertainty relation since A{tpt) = —ih/2 from 

C5±.) 

Let S-t be the surfaces defined by simultaneously imposing the constraints {Ci±} and 
the above gauge conditions. These are coordinatized by the physical variables — q, p, (Ag)^, 
A{qp), [ApY — and the one remaining gauge degree of freedom — t. It is straightforward to 
verify that the variables (At)^, A{tq) and A[tp) generate a Poisson ideal of the algebra of 
physical and gauge variables (i.e. the variables that do not involve pt). That is, {0j,X} 
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is a sum of gauge conditions with some coefficients. It follows that on S± the gauge-fixing 
conditions have a trivial Poisson algebra {(f)i,(f)j} ^ and a vanishing Poisson bracket with 
the remaining free variables. It is not difficult to see that Ci± remains ffist class on E±. 
Since E is a. function of the "physical" variables only, {E, 0j} pa and so 

{(f>i,Ci±}^{(l)i,Pt±E}^0. 

Furthermore, writing Csi = A{tpt) ± (p3p/E + const, one can quickly establish that Csi 
also remains first class but has a vanishing Poisson flow on E±. 

The remaining set of surface- defining conditions composed of C2±, Cs±, Ci± and 
{(j)i}i=i^2.'i is second class for all admissible values of the physical variables. This can 
be seen by relabeling the conditions as Xi = C'2±; X2 = Csi; Xs = C'4±! X4 = 0i) Xs = </'2) 
Xe = 03 and looking at the Poisson bracket matrix Ajj := {xj, Xj}- On S± the components 
of the matrix are 

/ ^{ih + 2A{qp)) ^{Apf 



Uh-A{qp) -(Ap) 



2 



-(Aqy -\ih-A{qp) 

-2ih 

^{ih + 2A{qp)) A{qp)-\ih {Aqf 

V ^{Apf {Apf \in + A{qp) 

Calculating the determinant one obtains the same result for both choices of the sign 

det[A] « -Ah' ^ + {A{qp)f + 2 [{Apf{Aqf - ^) {{Apf{Aqf - {A{qp)f) 

The determinant is non-zero in the region where reality, positivity and uncertainty condi- 
tions are imposed on the state — that is, if one demands 

(Ag)2,A(gp),(Ap)2eR 
(Ap)^ {Aqf > 

{Ap^iAqr - {A{qp)f > ^ . 

With these conditions in place, the sum of the terms inside the square bracket in the 
expression for the determinant is strictly positive, which means that the determinant itself 
is strictly negative. 

There is one important check that one needs to perform. The introduction of 0^ = 0, 
i — 1,2,3 makes the surfaces 'E± a mixture of first and second class and one is required to 
adjust the Poisson structure of the functions parameterizing the surfaces through the use 
of the Dirac bracket. Before we can identify q, p, (Aq)', A{qp), (Ap)^ as the expectation 
values and moments of a physical canonical pair of operators, we need to verify that 
their Dirac brackets on T,± are identical to the Poisson brackets one would obtain for the 
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quantum variables associated with a single canonical pair. The bracket may be computed 
as follows 

{/, ^?}D.ac ■■= {/, 9} - {/, X^} i^-'r {X„ 9} . (11) 

Using the fact that Aij (and hence also (A~^)*-') is off-block-diagonal and that the physical 
variables have vanishing Poisson brackets with the gauge fixing conditions, one can easily 
verify that the second term in the above definition vanishes for the brackets between q, p, 
{AqY, A{qp), {ApY (as well as t), and therefore their Poisson structure is unchanged as 
required by our interpretation. These variables are the remaining physical quantities up 
to second order on the reduced quantum phase space. 

To summarize, we impose the gauge-fixing conditions (pi = 0, i = 1,2,3, interpret q, 
p, {AqY, A{qp), {ApY as the physical expectation values and moments and as a result 
demand reality, positivity and quantum uncertainty. With all of these conditions taken 
together, 0j = 0, z = 1,2,3 restrict the gauge freedom up to the orbits generated by Ci± 
(recall that C5-1- generates no flow on This means that the time evolution is given by 

Cnam = fJ'lCl± ■ 

Finally, we fix the remaining Lagrange multiplier /zi by demanding t — the last remaining 
gauge variable — to be the time parameter. That is, we demand that {t, CHam} ~ 1, which 
leaves us with 

C-Ram = Ci± = pt± E . 

Taking the non-relativistic limit of E we recover the results for a deparameterized free 
Newtonian particle. Specifically, if we formally take p^ /ni^ to be of order 5, it follows that 
in a semiclassical state {ApY /m? is of order M. We expand the expression for E to the 
leading order in 5: 



m 



p^ + {ApfY / 4p2(Ap)^ 



E=^[l + ^- 1 + Wl- 



\ m? J \ y (j9^ + m? + {ApYY 

2m 

3.3 Comparison with the Klein-Gordon solution 

The standard positive frequency solutions to the Klein-Gordon equation (see for exam- 
ple [H]) form a Hilbert space of momentum-space wave-functions square integrable with 
respect to the Lorentz-invariant measure: 



7^ = lM R, — , where = \fW+ 
V 2efc / 
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The system can be described through a canonical pair of observables, represented on as 

■^(d fd 1 
p — k and q = «^ T^ + Cfc tttt; — 



^dk " \dk2e 

The time evolution is generated by the Hamiltonian H = (p^ + m^l)^. One can evaluate 
the evolution equations for the expectation values of the observables using Ehrenfest's 
theorem 

In our formalism, the right-hand side is equivalent to the quantum Poisson bracket between 
the expectation values, thus 

A(0) = {(0>, (H)} + ?m = {(O), ((p-^ + ,nH]i)} + ^. 

We recall that our procedure at order h together with the gauge fixing conditions for the 
positive frequency solutions resulted in 

^(O) = {(0),p, + E}^ {(O), E} + ^. 

In order to see whether the methods agree, we only need to compare ((p^ + m^l)^) and 
E to order h. To verify this explicitly we expand the operator in terms of its moments, 
assuming the expectation value to be taken in a semiclassical state: 

p^ + mH)^\ = ( (p^ + m^)^ H r(P - p) H — t(p - p)> + (higher moments) 



For comparison, we expand E in powers of (Ap)^, which we assume to be of order h. 



= V^^T^^i + J^'I'^J^, + o HAp)') 

Thus, up to the terms of order h the two results agree. 

Unlike the exact Klein-Gordon solution, our approach avoids explicit reference to a 
representation. The action of the Lorentz group on our variables can be understood through 
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its action on the algebra of observables. In particular, we assume that the pairs (pt, p) 
and (t, q) transform as components of a contravariant and a covariant vector respectively. 
Looking at the truncated system of constraints (ISl) under a Lorentz transformation one 
finds that C remains invariant, while the pairs (Cp„ Cp) and (Ct, Cq) themselves transform 
as components of a contravariant and a covariant vector respectively, so that the whole 
truncated system of constraints is preserved. 



3.4 Free massless particle 

For a massless particle, the constraint operator takes the form 

C = P?-P^ 

To second order in moments, the constraint functions produced remain as in equation ([H]), 
except for C, which now reads 

The disappearance of a constant term from C does not affect the Poisson algebra of the 
constraints, so that the table of Section 13.11 still applies. The solution to the constraints, 
however takes on a simpler form: following the same steps as previously and eliminating 
(Api)^ in a way compatible with the semiclassical approximation we obtain 

Together with C = this implies 

2 2 
Pt=P ■ 

As we see, the classical constraint is satisfied exactly by the expectation values. We solve 
this via 

Pt = ±\P\ ■ 

There are two related reasons for taking the absolute value of p in the above solution. 
Firstly, to emphasize the sign of the energy. Secondly, to match the limit as m is set to 
zero of the solution obtained in Section 13.11 The full solutions read 

Pt = ±\P\ 

(Ap.Y = (Apf 
A{ptq) = ±^(^A(gp) + y 

A{ptp) = ±^(AP)'- (12) 



13 



The steps of Section 13.21 can be repeated exactly for the m = case with \p\ playing the 
role of E. With the gauges fixed in an identical way, this results in evolution on q, p, 
(Ag)^, A{qp), (Ap)^ generated by the constraint 

Cnam =Pt± \p\ ■ (13) 
The implications will be discussed further in the conclusions. 

4 Relativistic particle in a potential 

In this section we consider the consequences of adding a potential term to the quantum 
constraint. We consider a quadratic time- independent potential in Section 14.11 followed 
by a homogeneous time-dependent potential in Section 14.21 The systems considered in 
this section have the same kinematical degrees of freedom as the free relativistic particle; 
however, the additional terms in the constraint element break Lorentz invariance. On the 
other hand, certain structural properties of the constraint element remain very similar to 
the free particle case, which makes extension of the calculations performed in Section [3] 
fairly simple. As we will see, the constraints are still straightforward to solve, but their 
Poisson algebra is only approximately closed and in the case of the time-dependent po- 
tential, the gauge analysis requires more subtlety. These examples show that the effective 
procedure used here is feasibly applicable to a wider range of models than the existing 
explicit constructions of a physical inner product. 

4.1 Quadratic potential 

A relativistic particle in a quadratic potential is subject to the constraint 

C 2 2 2 2-<i 

= Pt ~ P ~Q —ml. 

(A coupling constant in the potential could be absorbed by rescaling.) This gives rise to 
the following set of constraint functions truncated at second order 

C = p2 _ p2 _ ^2 _ ^2 _^ (^^^^2 _ (^^^2 _ (^^^2 ^ Q 

Ct = 2ptA{tpt) + ihpt - 2pA{tp) - 2qA{tq) = 
Cp, = 2pt{Aptf - 2pA{ptp) - 2qA{ptq) = 
Cg = 2ptA{ptq) - 2pA{qp) - ihp - 2q{Aqf = 

Cp = 2ptA{ptp) - 2p{Apf - 2qA{qp) + ihq = . (14) 

The above system of constraints is first class only to order h as can be seen from their 
Poisson algebra in Table [2l 

The system of constraints may be solved following the same steps that have been 
employed to solve the constraints for the free particle. We use C, Cg, and Cp to eliminate 
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Table 2: Poisson algebra of constraints for the particle in a quadratic potential. First terms 
in the bracket are labeled by rows, second terms are labeled by columns. 





c 










c 





— 2Cpj 







-2Cq 








4ptCpt-2pCp-2qCq 


2piC,+2pCt 


2ptCp+2gCt 


Ct 







+4A{ptp)A{tq) 


+4A(tg)(A(<7p)+f ) 


+4A(f<?)(Ap)2 








-4A{ptq)A{tp) 


+4(Ag)2A(tp) 


-4A(tp)(A(gp)-f ) 






-4ptCp( +2pCp+2qCq 




2pCpj 


-2qCp^ 


'-'Vt 





-4A(ptp)A(tqr) 





-4{Aq)^A{ptp) 


+4(Apf A(jitq) 






+4A{pig)A(ip) 




-A{ptq){A{qp)-%) 


+4A(ptp)(A(<irp)-f ) 






-2ptC'q-2pCt 


-2pCp, 




2ptCpt-4pCp-4(jCg 




— 2Cp 


-4A(tg)(A{gp)+f ) 


+4(Ag)2A(ptp) 





+4((Ag)2(Ap)2-^) 






-4A(tp)(Ag)2 


+4A(ptg)(A(gp)-f ) 




-4{A{qp)f 






-2ptCp+2qCt 


2qCpt 


4pCp+4gCg-2ptCpt 








-4A{tg)(Ap)2 


-4(ApfA(ptq) 


-4((Ag)2(Ap)2-/^) 









+4A(tp)(A(gp)-f ) 


+4A(ptp)(A(gp)-f ) 


-4(A(gp))2 





p^, A{ptq) and ^{ptp) respectively. Substituted into Cp^ this yields a quadratic equation 
in {Aptf: 

= {{AptfY - (p2 + g2 ^ ^2 ^ (Ap)2 + (AqY) {Aptf + {p\Apf + 2qpA{qp) + g'(Ag)2) 

The solution compatible with the semiclassical approximation has the form 

Pt = ±E 
A{tpt) = ^ (pA(ip) + gA(tg)) - ^ 

{Aptf = p2^^2^^2^(^^)2^(^^)2_^2 

±1 / ih 
A(ptp) = ^ (^p(Ap)2 + gA(gp) - ) (15) 

where 



E ^ — 
V2 



p^ + q^ + m^ + {Apf + {Aqf (16) 



+ (p' + + + (Ap)' + (A?)2)2 - 4 {p'iApf + 2?pA(gp) + ^'(Ag)^) 



Wc note that as C — pj — and pt are both exactly first class, E must also be first class, 
since it does not vanish on the constraint surface. 



15 



If one applies the semiclassical approximation once again to drop the terms of orders 
higher than h, the constraint system may be treated as first class. The linearized versions 
of the constraints take the form 



Ci± = pt±E 

C2± = {^pt? -p" -q^ -m" - {^pf -{^qf + 
C3± = A{ptp) ± i (^p{Apf + qA{qp) - gy^ 

1 / ill \ 

Ca± = A{ptq)±-[pA{qp)+p- + q{Aqf\ 

= A{tp,)±^{pA{tp) + qA{tq)) + '-^. (17) 

One can follow the process outlined in the Section 13.21 and impose the set of gauge-fixing 
conditions ffTOj) . Once again Ci± is first-class on the gauge- fixed surfaces S-t and writing 
C^± = A(tpt) ± §03 ± §02 it is not difficult to see that, once again, C5± remains first-class 
but has a vanishing flow on S-t. We recall that the gauge conditions have a vanishing flow 
on the remaining free variables; therefore only the first term in the expressions for each 
of the constraints Ci± above has a non- vanishing Poisson bracket with the conditions 0j. 
As a result, the Poisson bracket matrix A remains as in Section [3l2] up to entries of order 
h. Imposing reality, positivity and quantum uncertainty and demanding {t, Cnam} ~ 1 we 
once again obtain 

Cuam = Pt^E . 

Directly expanding E in powers of (Ag)^, A{qp) and (Ap)^ we get the expression to order 

h 



E = a/p^ + + 



^ (g2 + rn^){Apf - 2qpA{qp) + {p^ + m'^){AqY 



2(jP' + g^ + m?Y 
+ O {{Apr) + O ((Ag)^) + O ((A(gp))2) . (18) 



The specific constraint considered in this section can be implemented quite completely 
at the level of physical states, although specifics of the dynamics are more difficult to 
extract than with effective methods. Below we briefly describe the solution and compare 
it to our effective treatment. The algebra elements may be represented kinematically as 
differential operators on the space of square-integrable wave-functions in two variables xq 
and Xi in the usual way 

h d h d 



xo, Pt = -T^, q = xi, p = - 



i dxo i dxi 

The constraint operator splits into a sum of commuting, and therefore simultaneously 
diagonalizable, components: 
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• Pj = —^4^ has infinite-norm eigenstates of tlie form exp(^a;o)0(xi), with eigenval- 



• p^+q^ = —fi?-^+xl = 2Hharm5 where Hharm is precisely the standard Hamiltonian of 
the harmonic oscillator on Xi (with mass and frequency set to unity). This operator 
has normalizable eigenstates of the form ip{xo)ipn{xi), where ipn{xi) is the usual 
normalized n-th eigenstate of the harmonic oscillator, the corresponding eigenvalues 
are 2(n + ^)h. 

• Every wavefunction is an eigenstate of Im^ with the eigenvalue m^. 

The eigenfunctions of the constraint operator are '^k,n{xo,xi) = exp{^xo)'fn{xi), the 
corresponding eigenvalues are (/c^ — 2{n + — ni^). The space of solutions to the quan- 
tum constraint equation is therefore spanned by the wavefunctions '^k,n for which k = 

±'\j2{n+ + w?. These states are not normalizable with respect to the square integra- 
tion in both xq and xi, however they have unit norm when the integration is taken with 
respect to Xi alone. 

As is usually done for such systems we decompose the solution space into two segments 
one belongs to the positive part, the other one to the negative part of the spectrum of 
Pi (e.g. the separation of positive and negative frequencies of the solutions to the Klein- 
Gordon equation). The general element of the solution space is a linear combination of 
either positive or negative frequency null eigenfunctions of C, denoted by and \1'^ 
respectively 



The separation into two components allows us to define a positive-definite physical inner 
product on each one of them individually 



where the bar denotes a complex-conjugate and both states \E', $ belong to the same 
component. On the space of solutions, the above inner product is independent of the 
value taken by xq. Furthermore this inner product is consistent with interpreting t as 
time, since we can formally write (\l'|t|\E')phys = xq. (This equation requires some care in 
its interpretation since t is not a physical operator. But just viewing the integration on 
the right hand side of f|T9l) easily allows us to introduce an operator t by multiplication. 
The dependence on xq of the result is then in agreement with the fact that t is not a 
physical observable.) It is also consistent with the gauge choices of equation flTUl) . which is 
straightforward to verify using the fact that for any operator A, polynomial in t, pt, q, p 







ues k"^. 





(19) 



(^|tA|^)phys = Xo(^| A|^)phys . 
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The physical states may be interpreted as solutions of one of the two Scrodinger equations 



h d 



= ± (p2 + + mH) ' ^ 



That is, an ordinary quantum mechanical system with time evolution in the variable xq gen- 
erated by the self-adjoint, positive square-root Hamiltonian H = (p^ + q^ + m^l)^, which 



is defined through its action on the basis of eigenstates: Hy9„(a;i) = \/2(n + + m'^ipn{xi). 



To compare the physical states and the effective solutions, we proceed as we have done 
before, in Section [331 We expand the expectation value of the square-root hamiltonian in 
a semiclassical state 



Proceeding with the above expansion and keeping only the terms up to order h one obtains 
the expression that is identical to the one for E in equation f[T5]) . Thus, to leading order in 
the semiclassical regime, the effective solution to the constraint is consistent with physical 
state evolution, and the gauge choice of equation (fTOj) is consistent with the physical inner 
product defined above together with the interpretation of (t) as measuring the physical 
time. For a direct comparison between fully quantum and effective time evolutions for a 
specific semiclassical state of this system see App. [Bl 

If we replace (p^ + q^ + m^l) in the constraint by any positive operator /(q, p) analytic 
in q and p, physical states can in principle be found in a similar way. One could find the 
spectrum of /(q, p) and construct the solutions out of its simultaneous eigenfunctions with 
Pj . Finding the spectrum of a given operator is in general a complicated task. Further, we 
were helped in this example by the fact that the spectrum of (p^ + q^ + m^l) is discrete 
and the eigenfunctions are normalizable with respect to square integration over xi alone. 
Defining the physical inner product is more difficult if parts of the spectrum of /(q, p) 
are continuous. Finally, determining suitable coherent states for semiclassical purposes, as 
done for this model in App. [Bl can be a difficult task. The leading order effective solution, 
on the other hand, can be obtained in much the same way as was done for the above 
example, without explicit in-depth knowledge of the spectrum of q or the exact form of its 
eigenstates. 

4.2 Time-dependent potential 

As mentioned in the introduction, time- dependent potentials are of interest in quantum 
cosmology. Another example where time-dependent terms arise is a relativistic particle 
moving on a non-static curved background space-time. In such cases, our methods can be 





' + mH) - (p2 + q2 + mH) 
2(p2 + q2 + m2l)^ 



+ (higher moments) . 
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used as well, but additional subtleties do arise. Adding a "potential" V(t) = \t to the 
classical constraint gives the second order quantum constraints 

C = p^^~p^-m^ + (Apt)^ - {ApY + At = 
Ct = 2ptA{tpt) + ihpt - 2pA{tp) + \{Atf = 

Cp, = 2pt{Aptf - 2pA{ptp) + \A{ptt) - U\h = 
Cq = 2ptA{ptq) - 2pA{qp) - ihp + \A{qt) = 

Cp = 2ptA{ptp) - 2p{Apf + \A{tp) = . (20) 

These constraints, once again, form a closed Poisson algebra only up to order h. We 
proceed to solve the above set of polynomial equations explicitly — PtCp^ = implies 

1 ih 
Pti^PtY - PtP^iPtp) + -XptAitpt) - —Xpt = . 

Using pCp = and XCt = to eliminate ptpA{ptp) and XptA(tpt) respectively we obtain 

p',{Ap,Y - p\Apf + XpA{tp) - \x\Atr - ^Ap, = . 

Finally, we eliminate (Apt)"^ using C = to obtain a quartic equation in pt 

=pf - (p2 + rn''-Xt + (Apf) p] + '-^Xpt + {p\Apf + h.\Atf - XpA{tp)^ . 

The exact solutions to the above quartic equation are of course readily available, though 
they involve long algebraic expressions and are not particularly illuminating. Furthermore, 
due to the linear term in pt in the equation, the gauge choices we have employed previously 
together with the reality conditions would lead us to conclude that pt is complex. A more 
subtle gauge analysis is required to solve the constraint without further approximations. 
For instance, to respect reality conditions one would use a gauge relating moments to the 
expectation values, for example by making A{tp) dependent on pt. 

However, assuming the potential changes very slowly allows one to move forward with 
the standard gauge choice. Treating A as a second small parameter in addition to h}^"^, 
such that A^ is of higher than second order and discarding the terms of order higher than 
h we are left with 

pt - [p^ + - At + {Apf) pI + p^{Apf = . 

This equation could also be obtained by directly dropping products of A and second order 
moments or h in the expressions for the constraint functions (!20|) and solving them. It is a 
quadratic equation in pf, with solutions that are much easier to interpret. Compatibility 
with the semiclassical approximation once again selects for us a set of solutions that have 
a very similar form to those in ([7]), (IT^ . and (IT^ . with 



V2 



p'^ + -Xt + {Apf + + - At + {Apf f - Ap'^^Apf) ' 
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One can then repeat the gauge-fixing procedure we have previously employed and recover 
Cnam = Pt ^ E, for a slowly varying potential in a semiclassical state. Thus, to the 
semiclassical order considered, the system behaves as a non-relativistic quantum particle 
in one dimension subject to a time-dependent hamiltonian H = (p^ + m^l — Xt) ^ . 

In fact, a more general "slowly varying" potential may be treated to order h in an 
analogous manner. We assume that the potential has the form V{t) = V{0)1 + XV{t), 
where A is small in the sense discussed earlier, and V{t) is a polynomial in t. This implies, 
in particular, that (AV^(t)) = \V(t) + 0{h2). We absorb the constant part of the potential 
into the definition of m and the constraint functions to order h look exactly as they do for 
a free relativistic particle ([5]), except for C, which acquires an extra term 

The resulting system of constraint functions may be solved, gauge-fixed and interpreted 
directly following the methods employed throughout Sections |3] and HI This demonstrates 
the fiexibility of the constructions, confirming the methods of pTl[T2l[T3] . where also slowly 
varying potentials were assumed. In contrast to this earlier work, the general methods 
presented here are in principle applicable to arbitrary potentials, but the gauge fixing 
would have to be considered in each case in detail. This provides access to questions about 
the role of time when potentials forbid a global monotonic internal time t, resulting in a 
new perspective to be followed elsewhere. 

5 Conclusion 

One of the main hurdles for quantum gravity is the physical Hilbert space issue. At least 
for semiclassical questions, technical and conceptual difficulties can be circumvented by 
using expectation values and moments directly rather than states. Other advantages of 
this method are that the specification of semiclassical regimes is easier via moments (while 
semiclassical wave functions are often difficult to formulate, even simple-looking Gaussian 
ones not always being semiclassical at all in some models of quantum cosmology as pointed 
out in [IS]) and that density states are automatically included. 

We have extended the methods for effective constraints of [2] to relativistic systems, 
clarifying several physically relevant issues of effective equations: 

• Square root effective Hamiltonians, which so far were strictly justified only for time- 
independent potentials, are valid even in the time dependent case provided the po- 
tential varies slowly in time. No extra conditions on the dynamics are implied by the 
positivity conditions. 

• Quantum back-reaction follows reliably from square-root Hamiltonians. Our exam- 
ples of relativistic systems have shown three different cases: 

— Massless particles do not contain moments in their reduced Hamiltonian flTS]) 
and thus are not subject to quantum back-reaction. 



20 



— Free massive particles do have quantum back-reaction from moments of all or- 
ders, which is initially unexpected since the quantum constraint has only a linear 
term of (Ap)^ with no coupling to the expectation values. (Our formulas, done 
here only to second order, do not show this explicitly.) 

— Particles subject to a g-dependent potential receive quantum back-reaction from 
the covariance of their wave function as seen in ( fT6l) . This is also unexpected 
since the expectation value of the constraint does not contain mixed moments. 
As in the case of massive particles, the unexpected results are explained by the 
presence of higher order constraints. 

• Higher order constraints, which are crucial for the effective procedure, also affect the 
amount of spreading of wave functions, or the time dependence of moments. For a 
free, massless particle wave packets do not spread, but they do in the other cases. 

Effective equations obtained in the way developed here reliably describe the physical be- 
havior of dynamical wave packets. With these considerations, effective methods as they 
have been used in quantum cosmology are established even in the case of 0-dependent 
potentials, as studied specifically for instance in [I3l [8] . A further application would be 
to reconsider the appearance of certain future singularities which have been shown not 
to be removed by the tree-level approximation (disregarding all moments) [20] but where 
quantum back-reaction is expected to be strong. 

Especially in the presence of potentials, calculations performed here are much more 
feasible than the methods involving constructions of physical Hilbert spaces followed by 
computations of the expectation values in explicit physical states. They can be expected to 
be of far more general applicability, including full quantum gravity. For such an extension, 
several other issues remain open: describing field theories and handling situations of many 
classical degrees of freedom. But there is already a promise that effective techniques allow 
one to evade difficult obstacles from physical inner product issues which so far have impeded 
progress. Especially the semiclassical regime of canonical quantum gravity and potentially 
observable effects can be brought under much higher control. 
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A Second order Poisson algebra 

The expectation values obey the classical Poisson algebra for two canonical pairs, where 
the non trivial brackets are 

{t,pt} = l and {q,p} = l. 
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Table 3: Poisson algebra of second order moments. First terms in the bracket are labeled 
by rows, second terms are labeled by columns. 





{Atf 


A(tpt) 


(Apt)' 


(A<7)2 


A(gp) 


(Ap)2 


A{tq) 


A(ptp) 


A(tp) 


A(ptg) 


(At)2 





2(At)2 


4A(tpt) 














2A(tp) 





2A{tq) 


A(tpt) 


-2(At)2 





2(Apt)2 











-A{tq) 


A(ptp) 


-A(ip) 


Mptq) 


(Apt)' 


-4A{tpt) 


-2(Api)2 














-2A{ptq) 





-2A(ptp) 





(Ag)2 














2(Aq)2 


4A(qp) 





2A(ptq) 


2A(tg) 





A(qp) 











-2(Ag)2 





2(Ap)2 


-A(tg) 


A(ptp) 


A(tp) 


-A(ptg) 


{Apf 











-4A(<?p) 


-2(Ap)2 





-2A(tp) 








-2A{ptp) 


Aitg] 





A(ig) 


2A(ptg) 





A{tq) 


2A{tp) 





A(tpt) 
+A(gp) 


(At)' 


(Ag)' 




-2A(tp) 


-A(pip) 





-2A(ptg) 







-A(tpt) 
-A(gp) 





-(Ap)2 


-(Apt)2 


A(tp) 





A(tp) 


2A(ptp) 


-2A(tg) 


~A{tp) 





-(At)' 


(Ap)' 





A(gp) 
-A(tpt) 




-2A(t(j) 


-A(ptg) 










2A(ptp) 


-{Aqf 


(Apt)2 


A{tpt) 
-A(<jp) 






The Poisson brackets between the expectation values and the moments are zero. In Ta- 
ble [3] we provide the Poisson brackets for the second order moments of quantum variables 
associated with two canonical pairs t, p^; q, p. 

B Coherent state and effective evolution 

Here we consider the constrained system of Section 14.11 in a specific semiclassical state and 
compare the evolution of "classical" quantities given by the effective quantum theory to 
properties of the state. This appendix is included to address the riliability of an effective 
solution to a constrained quantum system through comparison within a specific example. 
Recall the constraint: 

C = p2 - p2 - q2 _ mH . 

Formally, a positive frequency solution of the constraint described in Section 14.11 reduces 
the system to one canonical degree of freedom that evolves subject to the Hamiltonian 
H = (p2 + + m^l)^. We begin by providing the classical trajectory: a canonical pair 

q, p subject to the Hamiltonian function H = (p^ + + m^)^ evolves according to the 
equations of motion: 

P {p^ + q'^ + m^) ' = pH^^ 
-q (p^ + + m^)~' = qH~^ . 
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Using the fact that H itself is a constant of motion we differentiate the first equation above 
with respect to time to obtain the second order equation 

This equation has the general solution 

q{t)=Asm (^] +5cos 



where A and B are constants. It follows that 



pit) = H^q = Acos(^\ - sin 



and thus H = y/ A'^ + B'^ + iv?. The classical phase-space trajectory is therefore a circle of 
radius \/ A^ + B"^ traversed with the angular frequency {A^ + i?^ + m^) ^ . 

Effective equations of motion are set up in much the same way as their classical coun- 
terparts. Phase-space degrees of freedom are g, p, (Ag)^, (Ap)^, A(gp). The time evolution 
is generated by the function E of equation (ITS!) through the quantum Poisson bracket 

d _ p p(Ag)2 (2g2 - p2 - m^) + gA(gp) (4p2 - 2g2 - 2m2) - 3p(Ap)2 (g2 + m^) 

dt^ + g2 + ^2 2(p2 + g2 + ^2)f 

d — g 3g(Ag)2 (p^ _|_ ,77,2) _ pA{qp) (4g2 — 2p^ — 2m^) — g(Ap)2 (2p2 — q^ — m^) 

d/ ^p2 + g2 + ^2 2(p2 + g2 + ^2)f 

d ^^^^2 _ 2A(gp) (g2 + ^2) - 2(Ag)2gp ^^^^^2 _ 2(Ap)2gp - 2A(gp) (g2 + ^2) 



dt (^p2 + g2 _^ ^2^ 2 dt (^p2 + g2 _^ ^2") 2 

d ^ , , (Ap)2 (g2 + ra^\ _ (Ag)2 (p2 _^ 7„2n 
— A(gp) = 3 . 

dt (p2 + g2 + ^2)2 

The first term in each of the first two equations is identical to the classical equations of 
motion, the extra terms constitute the leading order quantum corrections. The system 
of equations is straightforward to evolve numerically for a sufficiently short time starting 
from a state that is initially semiclassical. 

The quantum evolution of positive frequency solutions is governed by the Schrodinger 
equation 

^ d i 

-^{x, t) = (p2 + + m^l) 2 ^(x, t) . 

i at ^ 

The square root operator has the same eigenstates {(Pn{x)}n=o,...oo; as the one dimensional 



quantum harmonic oscillator, with eigenvalues = y2{n+ ^)h + m?. A wavefunction 
can be evolved by decomposing it into the eigenstates of the Hamiltonian. If the state at 
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i = is given by ^(x, 0) = Yl'^=o^nfn{x), where c„ are constant complex numbers, then 
at any other time, the wavefunction is 



= ^c„exp (^-^-^^ Vn{x) . 



(21) 



To compute the expectation values we write q = ^Jhj^ia* +a), p = i\Jh/2 (a* — a), 
where d* and d are the usual creation and annihilation operators of the quantum harmonic 
oscillator. One finds 

c„c„+iexp -it- 



n=0 
oo 

(*,p*)(0 V2%+l)Im 



n=0 



CnCn+1 exp -it 



h 
h 



In a similar manner one can obtain expressions for the moments of q and p. 




Figure 1: Classical (dotted), coherent state (solid) and effective (dashed) phase space 
trajectories, evolved for < i < 5go- 

In order to complete the comparison we select a semiclassical state with a known de- 
composition into the cigenstates {ipn{x)}- A simple choice is to set the initial wavefunction 
to be a coherent state of the harmonic oscillator 



Cn = exp 



\a\ \ a 



(22) 



which we can consider as a kinematical coherent state for our system. For the non- 
relativistic harmonic oscillator, as time goes on a changes, but the shape of the state 
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is preserved. Clearly, this is not the case for our relativistic evolution: Combining (!22|) 
with (12T|) . we have an evolving state which, expanded in harmonic oscillator stationary 
states, has coefficients 

Due to the square root in this relativistic model, these coefficients are not of the coherent 
state form (|22|) unless t 7^ 0. The physical states we obtain are not dynamical coherent 
states; quantum back-reaction ensues which in the effective treatment is captured by the 
coupling terms between moments and expectation values in f|T8|) . 

For a specific numerical example, we set a = so that at t = the state is a Gaussian 
peaked about q = qo and p = 0, with zero covariance and minimal spread 

(Agf = (Ap)2 = ^, A(gp)=0. 

We take these as initial values for the numerical evolution of the effective equations for the 
system. We assume the two physical scales to be separated by a single order of magnitude 
by setting = 10. For simplicity we set m = 0. Depicted in FIG. [1] are the classical, 
coherent and effective quantum phase-space trajectories starting from the same initial 
state. One can see that the effective equations describe the correct semiclassical trajectory 
for much of the evolution displayed. An internal measure of consistency is the size of 
second order moments. From FIG. [2], we see that the semiclassical approximation clearly 
breaks down after approximately t = 2qo, as (Ap)^ becomes too large. The same figure 
demonstrates that up until that time the evolution of the moments themselves is very well 
approximated by the effective equations. The other two moments display similar behavior. 

Ap 



90 




t 



10 2qo 3qo 4qo 



Figure 2: Coherent state (solid) and effective (dashed) evolution of the second order mo- 
ment Ap = a/ (Ap)2 in units of go- 
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